Edge-localized states in quantum one-dimensional lattices 
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In one-dimensional quantum lattice models with open boundaries, we find and study localization 
at the lattice edge. We show that edge-localized eigenstates can be found in both bosonic and 
fermionic systems, specifically, in the Bose-Hubbard model with on-site interactions and in the 
spinless fermion model with nearest-neighbor interactions. We characterize the localization through 
spectral studies via numerical diagonalization and perturbation theory, through considerations of 
the eigenfunctions, and through the study of explicit time evolution. We concentrate on few-particle 
systems, showing how more complicated edge states appear as the number of particles is increased. 

PACS numbers: 63.20.Pw, 63.20.Ry, 63.22.-fm, 03.65.Ge 



I. INTRODUCTION 

Experimental techniques in the fields of mesoscopics 
and ultracold atoms have advanced to the point where 
it is feasible to explore the physics of few-boson and 
few-fermion systems, in particular in lattice systems 
which have traditionally been the basis for many-particle 
physics studies. The study of a few quantum particles 
in lattice situations provides surprises and unexpected 
phenomena quite distinct from issues in bulk condensed 
matter physics, whose focus is on man^z-particle ground 
states, and from atomic or optical physics, where lattice 
systems are not very common. For example, a recent ex- 
periment has explored repulsive binding of boson pairs 
in a one-dimensional optical lattice [H. This is there- 
fore the appropriate time to investigate further intricate 
and non-intuitive phenomena involving a few quantum 
particles in lattice systems. In this article, we consider 
interacting bosons or fermions in one-dimensional finite 
lattices, and present studies of localization at the lattice 
edge. We characterize the edge-localized states, and dis- 
tinguish them from bound non-localized states, through 
analyses of the energy spectrum and band structure, den- 
sity profiles of eigenstates, and dynamics. 

For bosons, we use the well-known Bose-Hubbard 
model which has attracted a great deal of atten- 
tion in the last decade due to its relevance to describ- 
ing laser-cooled bosonic atoms subjected to an optical- 
lattice potential [1, 0] . Localization in the Bose-Hubbard 
model is of particular interest because its large-boson 
limit can for many purposes be approximated by the 
discrete nonlinear Schrodinger (DNLS) equation, which 
displays a host of localization phenomena. Nonlinear- 
ity allows time-periodic and spatially-localized solutions 
of the DNLS (and other lattice nonlinear differential 
equations) known as intrinsic localized modes or discrete 
breathers [^, 0, B @- More relevant to the present 
work, the DNLS on finite lattices also possesses edge- 
localized modes, sometimes called discrete surface soli- 
tons [l^, [III, [l^l . It is therefore expected that the large- 
boson limit of the open-boundary Bose-Hubbard model 
will possess eigenstates in which the bosons are localized 
at the edge. In this Article, we pose the edge-localization 



question in the extreme opposite limit of a Jew quantum 
particles, where the mean-field approximation (DNLS 
equation) cannot a priori be expected to provide the cor- 
rect intuition. The answer turns out to be subtle — this 
phenomenon is not present for the case of two particles, 
but appears when the particle number is three or more, 
as follows from numerical studies in Ref. fl3|. 

Remarkably, we find that the phenomenon is not re- 
stricted to bosons, but also happens in other quantum 
lattice models. For example, the spinless fermion model 
with nearest-neighbor interactions, sometimes known as 
the t-V model, has similar behavior, i.e., it possesses no 
edge states with two fermions, but does have such local- 
ized eigenstates with three fermions. Since this model 
does not allow multiple occupancies, localization in this 
case refers to a sequence of neighboring sites being oc- 
cupied, rather than all the particles clustered in one site 
as in the case of bosons. Although the current work will 
focus on the two models mentioned above, our finding in- 
dicates that edge-localization may well be a generic phe- 
nomenon in quantum finite lattice systems. 

Localization in quantum models generally does not ap- 
pear in simple Hamiltonians, but instead requires disor- 
der [14] or impurities breaking the translation invariance. 
(This is in contrast to lattice differential equations, e.g., 
the DNLS equation, where the interplay of nonlinearity 
and spatial discreteness is sufficient to create localiza- 
tion.) It is therefore of significant theoretical interest to 
explore the simple mechanism for quantum localization 
that is studied here, requiring only open boundary condi- 
tions in an interacting lattice Hamiltonian. We also note 
that localization due to impurities/disorder is generally 
a single-particle effect, while the mechanism we present 
is a collective phenomenon since it requires at least three 
particles. 

While we focus on small numbers of bosons or fermions, 
the basic message is that edge states exist for any num- 
ber of particles larger than two. We give an expla- 
nation based on perturbation theory. We also provide 
some hints toward the many-particle situation, by show- 
ing that the four-boson case has an additional type of 
edge-localized eigenstate in addition to the obvious gen- 
eralization. 
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FIG. 1: (Color online.) Energy spectrum of 10-site Bose- 
Hubbard chain with 7 — 10. Left panels (a,c) show the case of 
two bosons (n = 2), right panels (b,d) show n = 3. Top panels 
(a,b) plot ordered energies against index v, for fixed hopping 
parameter t = 1. Lower panels (c,d) plot energies against hop- 
ping strength t. Insets focus on the lowest (breather) band, 
showing that two edge-localized states separate out for three 
bosons (b,d) but not for two bosons (a,c). 

Model Hamiltonians. We consider one-dimensional lat- 
tices with L sites subject to open boundary conditions. 
The Bose-Hubbard model Hamiltonian is 

L L-l 

(1) 

Here a and are second-quantized bosonic operators. 
The first and second term describe respectively on-site 
attractive (7 > 0) interactions and nearest-neighbor hop- 
ping. The edge-localization physics is almost unchanged 
in the case of repulsive interactions. 

The model for spinless fermions is described by 

L-l L-l 

(2) 

where V is the (repulsive) nearest-neighbor interaction 
strength, and cj and cj are second-quantized fermionic 
operators. 

II. BOSE-HUBBARD CHAIN WITH TWO AND 
THREE BOSONS 

We will first describe the cases of two bosons and of 
three bosons in the Bose-Hubbard chain, and show how 
the two situations differ by virtue of the latter having 




FIG. 2: (Color online.) Spatial profile of site occupancies for 
several eigenstates of the three-boson Bose-Hubbard chain. 
Inset shows the same plot in semilog scale, the linear behavior 
indicating exponential localization in the edge states. Here 
L = 10 and 7 = 10. 



edge states. Some of the physics described in this section 
appears in Ref. [l3|. 

Strong interactions. — We first focus on large values 
of 7, which is the most relevant parameter regime for 
localization physics. 

The Hilbert space size (number of basis states) for n 

bosons m L > n sites is d„ = ^ ^ 'n ' ™ particular 

d2 = L{L + l)/2 and = L{L + 1)(L + 2)/6. Figure 
[T] displays spectral properties and band structure, i.e., 
the distribution of the d„ eigenenergies. We label the 
eigenstates from low to high energies with the label v 
running from 1 to d„. 

The most prominent feature of the large-7 spectrum 
is the band structure. For the n-boson system there are 
p{n) bands, where p(n) is the number of integer partitions 
of 71. In the two-boson case (Figure [T]-a and[l]-c), there 
are two bands. In the upper band around zero energy, the 
dominant contributions come from configurations where 
the bosons do not share the same site. The lower band 
contains the L two-boson bound states, dominated by 
linear combinations of configurations where the bosons 
sit on the same site. These are the quantum analogs of 
classical discrete breather solutions. This band is thus 
called the quantum breather band or the soliton band 
p^ . In these states, the separation probability of the 
bosons decays exponentially with distance [13, [iM El] ■ 

In the three-boson case (Figures[T]-b and[T]-d), the spec- 
trum contains three energy bands. The lowest-energy 
band is formed by the L three-boson bound states (three- 
boson breather band), where there is a high probability 
of finding the three bosons at the same lattice site. The 
second band from below is formed by L{L — 1) "2-1-1- 
boson states" , where there is a high probability of find- 
ing two bosons at the same lattice site with the third 
boson elsewhere. Finally, the third band is the three- 
boson continuum, whose L{L — \){L — 2)/6 eigenstates 
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FIG. 3: (Color online.) Size-dependence of max{nj) for (a) 
all eigenstates of the two-boson Bose-Hubbard chain, and (b) 
several eigenstates of the three-boson Bose-Hubbard chain la- 
beled by i/. Here 7 — 10, and L goes from L = 5 to L = 17. 



are characterized by having the three bosons in different 
sites. 

The edge states can be identified by zooming onto the 
three-boson bound state band (insets of Figures [T]-b and 
HJ-d). We note that two states stand out from the rest 
of the band, with larger splitting. These are the edge 
states, as we demonstrate further below. Because of re- 
flection symmetry, the dominant contributors to the two 
eigenstates are not the left-edge and right-edge states 
{\Ei) = |3000...) and \Er) = |...0003)) directly. Instead, 
the eigenstates are dominated by the linear combinations 

± \Er )^ . The remaining {L — 2) eigenstates of 

the breather band are dominated by linear combinations 
of the remaining {L — 2) three-boson bound state configu- 
rations, 103000...), 100300...), ... I...0030). These (L-2) 
states do not mix with the two edge states because of the 
energy splitting. 

Note that the energy-separated edge-localized states 
are not present in the two-boson case (insets of Figures 
dJa andHJc). 

To see the localized nature of the edge states, in Fig- 
ure [2] we plot density profiles (site occupancies), i.e., ex- 
pectation values of boson number at each lattice site, 
(nj)i. = {Xiy\a]aj\xu) , \xiy) being an eigenstate. Figure 



[5] shows 



against j for the edge states, and also for 



two other eigenstates (non-edge breather states) for com- 
parison. We see that the edge states are exponentially 
localized at the edges of the lattice, as seen through the 
linear behavior of In(nj) in the inset. 

The localization is further demonstrated through the 
scaling of the largest site occupancies, max(nj), with 
system size L (Figure [3]). For eigenstates that are lo- 
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FIG. 4: (Color onhne.) Spectrum and density profiles at 
smaller interactions, for three-boson system in 10 sites {t = 1). 
(a) The bands merge at small enough 7. (b) For this system 
the breather band loses its identity around 7 ~ 1.8. We follow 
the edge-localized states at smaller 7; shown surrounded by 
a thick-lined green polygon. At smaller 7, the v — (i — 1) 
and V = L states are no longer the edge states, (c) Site occu- 
pancies at 7 = 1.4. The almost exponentially localized curves 
are for the edge states, and the two non-localized curves are 
the 1/ = [L — 1) and v = L states, which are no longer the 
edge-localized states, (d) Size-dependence of largest site oc- 
cupancy. 



calized, this quantity will be independent of the system 
size, whereas for extended states it should be a linear 
function of 1/i because the bosons are spread across L 
sites in extended states. In the two-bosons case (Figure 
[3]-a) all eigenstates are extended, thus max(nj) depends 
linearly on 1 jL for all of them. In Figure [3]-b for three 
bosons, most states also have max(nj) varying as 
except for two states for which max(nj) are independent 
of L. This flat set of points in Figure [3]-b clearly illus- 
trates the localization phenomenon. 

We have described the spectrum for the attractive 
Bose-Hubbard model. The spectrum for the repulsive 
case (7 < 0) is obtained simply by inverting the energies 
(Ev^ — Ey). The eigenfunction characteristics described 
in Figures [2] and [3] are identical in the repulsive case. 

Weak interactions. — It is natural to ask whether 
the edge-localization survives for weaker interactions. In 
particular, for three bosons and 7 < 2, the breather band 
merges with the 2-1-1 band (Figure |3]-a). One might 
speculate that there might be a sharply defined value of 
7 below which there is no localization at the edges. 

We address this question by following the two edge 
states adiabatically to lower values of 7, for three bosons 
in 10 sites, as shown in Figure [H-b. Note that, after the 
breather band has merged, the v = L—l and v — L states 
are no longer the edge states. In Figure [3]-c we plot the 



spatial density profiles of the two edge states, as well as 
the non-localized = L—1 and v — L states. Figure |3]-d 
shows the size dependence of niax(nj). The localization 
is weaker than in the large-7 situation (niax(nj) is less 
than 1.5 and the exponential decay is imperfect), but it 
is still present. Thus the localization phenomenon does 
not completely disappear at some sharp value of 7, at 
least for small {L ^ 10) lattice sizes. 
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III. PERTURBATION THEORY 

The basic explanation for the edge-localization phe- 
nomenon is that the edge states split off from the rest of 
the breather band, for three or more particles. We will 
now explain this energy splitting through perturbative 
calculations in the hopping parameter. It is helpful to 
introduce an effective single-particle model, which con- 
tains only the parts of the Hilbert space relevant for this 
analysis. 

For the two-boson case, we keep the L breather band 
states, each corresponding to both bosons in the same 
site. We keep only those states of the continuum that 
are necessary as intermediate configurations to go from 
one breather state to another, i.e., states where the two 
bosons are in neighboring sites. In Figure [5]-a, the 2- 
boson states are shown with filled circles and the in- 
termediate l-|-l-boson states as open circles; together, 
they form the single-particle effective tight-binding chain. 
This effective chain has two different on-site energies 
£1 ~ —7 and £2 = alternating along the chain (Fig- 
ure [S]-a). The effective hopping strength is v^i, since 
the hopping is always from or to a doubly occuppied site 
in the original model. For a fixed size L of the Bose- 
Hubbard chain, the size of the effective single-particle 
chain, i.e., the number of states retained from the origi- 
nal Hilbert space, is L + (L — 1) = (2L — 1). 

Similarly, for the three-boson case we retain the L 
breather band states, and only those states of the (2-1-1) 
band that are necessary as intermediate configurations to 
go from one breather state to another, i.e., states where 
the two bosons and the lone boson are in neighboring 
sites. The single-particle chain now has two on-site ener- 
gies £1 = —37 and £2 = —7, and two hoppings ti = \/?>t 
and ^2 = They form a chain of basis three per unit 
ceh (Figure [Sl-b). 

Degeneracy splitting — At zero hopping (t = 0) or in- 
finite interaction, all bands, including the breather band 
which we are particularly interested in, are perfectly de- 
generate. We analyze the splitting of the breather band 
spectrum perturbatively in the hopping parameter t. The 
Hamiltonian \s H = Hq + tHi , where Hq is the interac- 
tion term. Since the relevant states of the Hilbert space 
are arranged along a single chain in the single-particle 
effective picture, we can use the position on the effective 
chain as a label for the relevant states. For example, for 
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FIG. 5: (Color online.) Construction of the effective single- 
particle chain for the Bose-Hubbard model with (a) two 
bosons and (b) three bosons. The sequence of hops of the 
two bosons (grey circles), initially at the same lattice site, are 
shown with the corresponding values of the matrix element of 
the Bose-Hubbard Hamiltonian. The energies ei and £2 of the 
states after each hop are also shown. The resulting effective 
single-particle chain is shown in the lower part of the figures. 



the two-boson system, 

L L-l 

Ho^J2 ei|2m- l)(2m- 1| + ^ £2|2m)(2m| , (3) 



while the perturbation Hi contains hopping terms like 
\2m){2m+ 1|. 

For two bosons, the lowest order at which the pertur- 
bation has nontrivial effects is 0{t^). The hopping per- 
turbation at this order already connects the states of the 
ground state manifold to each other, leading to a com- 
plete lifting of the degeneracy, so that the non-degenerate 
breather band emerges. The split energies are approxi- 
mately E2m-1 « -7 - 4^2(1 -I- COSkm)/l- 

For the three-boson case, one can carry out a similar 
analysis. The crucial difference is that, at second order 
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FIG. 6; (Color online.) Energies from second-order pertur- 
bation theory (dashed and dotted lines) compared with full 
spectrum (solid lines), for 7 = 10 and L = 10. 

in t, the lowest-manifold states are each connected to 
themselves, and thus receive an energy shift, but are not 
connected to other states within the manifold, and there- 
fore the degeneracy is not lifted. The energy shifts are 
different for the edge and non-edge states, because there 
is a single path for an edge state to couple to itself via two 
hopping events, while each non-edge state has two such 
paths. This is visually obvious through hopping events 
in real space in our effective single-particle chain (Figure 
[5]). The shifts at second order are 

£^cdgc = —37 — —t^ , 

-^non — edge ^ 3'^' t • 

7 

At the next order [t^) the degeneracy of the non-edge 
states is lifted, since three hopping events are required to 
connect two distinct members of the breather manifold. 
The degeneracy of the two edge states is only broken at 
much higher order. 

This analysis reveals the reason for the separating out 
of the edge states from the rest of the breather band. The 
energy shift of the edge states happens at lower order in 
the hopping than the order at which the degeneracy of 
the breather band is lifted. Therefore the edge states are 
robustly separated out for large interactions. 

Figure [5] plots the second-order perturbative results 
for the breather band energies, comparing with the ex- 
act energies computed numerically. The splitting of the 
edge-localized states from the rest of the band, which is 
the essential issue here, is well described by perturbation 
theory. The degeneracy lifting of the non-edge breather 
states is not captured in the second-order expressions. 

IV. FOUR OR MORE BOSONS 

Edge-localized eigenstates also exist in Bose-Hubbard 
chains with n>3 bosons; the perturbative argument we 
provided for 3 bosons can readily be extended to the 
general case. At strong interactions, the lowest band is 
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FIG. 7: (Color online.) (a) Energy spectrum of the Bose- 
Hubbard chain with four bosons, plotted against eigenvalue 
label. The inset zooms onto the breather band, while (b) 
zooms onto the (2-|-2)-boson band, as indicated by the dotted 
red arrows. Here 7 = 10, t = 1, and L = 10. (c) Spectrum 
of same chain as a function of the hopping parameter t. The 
insets focus, from left to right, on the 4-boson bound state 
band, the 2+2-boson band, and the 2-|-2-boson bound state 
sub-band, showing the separation of edge-localized states out 
of the former and the latter. 

the breather band with n-boson bound states. At t = 0, 
this band is collapsed as the L-fold degenerate ground 
state. For i ^ the degeneracy is lifted at order ©(i"), 
but the edge states already have a distinct energy shift 
at 0{t^), leading to edge localization. As in the n = 3 
case, one can also visualize the perturbative calculation 
with an effective single-particle model retaining only the 
relevant basis states. This will now have rt-site unit cells, 
i.e., sites representing breather states separated by 71 — 1 
sites representing states from other bands. 

In addition to the edge-localized states with all n 
bosons situated at the edge, for n > 3 the open Bose- 
Hubbard chain also has edge states with more compli- 
cated structure. We demonstrate this for n = 4 in Fig- 
ure [71 Other than the two edge modes on top of the 
breather band, we see features in the 2+2 band. As in the 
translation-symmetric case [20| , this band has a sub-band 
separating out (middle inset in Figure[3-c). This subband 
of L — 1 states is characterized by two doubly-occuppied 
sites neighboring each other {2-\-2 bound states), while 
the rest of the 2-1-2 band is dominated by two occup- 
pied sites at larger distances from each other. Unlike 
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FIG. 8: Profile of tfie site occupancy in tlie Bose-Hubbard 
cliain with four bosons for four eigenstates in (a) tfie 4-boson 
bound state band and (b) tlie 2+2-boson bound state band. 
Tfie same plots in semilog scale are shown in the insets, show- 
ing exponential localization. Here L = 10 and 7 = 10. 



the translation-symmetric case, however, if one zooms 
in further onto this sub-band, two edge-locahzed states 
separate out (rightmost inset in Figure [7]-c) . These edge 
states have the structure of two bosons at the edge, and 
the other two at the next-to-edge site. Figure [5]-a and 
[5]-b show the density profiles of these two types of edge 
states in the n = 4 case. 

The above results demonstrate the existence of more 
and more complicated additional edge states as the num- 
ber of particles is increased. 



V. SPINLESS FERMIONS 

We now turn to the spinless fermion model described 
by the Hamiltonian Figure [H displays through 

numerically calculated spectral properties that in this 
model, edge states do not exist for the two-fermion case, 
but appear when there are three or more fermions. The 
situation is thus similar to the Bose-Hubbard model. 

We first note that the spectrum of this model con- 
tains a breather band as in the Bose-Hubbard model [lB| ■ 
Fermi statistics forbids multiple occupancy of the sites; 
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FIG. 9: (Color online.) Energy spectrum of the spinless 
fermion model, for two fermions in the left panels (a,c) and 
three fermions in the right panels (b,d). Here L = 15 and 
V = 10. Insets focus on the (topmost) breather band. The 
three-fermion cases (b,d) have two states separating from this 
band; the two-fermion cases (a,c) do not have this feature. 

SO the breather modes correspond to all the fermions clus- 
tered in a connected segment of the lattice. Since we are 
using repulsive interactions {V > 0), the breather band 
now appears at the top of the spectrum. 

The spectral splitting of edge states is analogous to 
what we have described in the Bose-Hubbard case. For 
the two-fermion case, the breather band is completely 
split because the separation for edge and non-edge states 
all occurs at the same (second) order in t/V. For three 
or more fermions, however, the two edge states split off 
from the main band at lower order than the degeneracy- 
lifting of the rest of the band. As a result, there are now 
two edge-localized states separated at the bottom of the 
breather band. 

Figure [TO] shows density profiles for the edge- localized 
states as well as some non-localized states from the 
breather band. Each edge state now has a 'width' of 
three sites, because of fermionic statistics forbidding dou- 
ble occupancies. Since the two eigenstates are predom- 
inantly linear combinations of left-edge jlllOOOO...) and 
right-edge |... 0000111) states, they have occupancies of 
njWl/2 at the last three sites of each edge. As in the 
Bose-Hubbard case, the logarithmic plot makes clear the 
exponential nature of the localization. In this case the 
exponential decay starts after the third site. 

In Figure[Tll we display the robustness of the 3-fermion 
edge states by plotting the size-dependence of the max- 
imum site occupancy, for a fixed interaction potential. 
Analogous to the Bose-Hubbard model (Figure [S]) , the 
two-fermion case has max{{nj)^} values all varying as 
~L~^, for all states. In the three-fermion case, all 
eigenstates except two also have linear 1/L-dependence, 
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FIG. 10: (Color online.) Spatial profile ol fermion occupan- 
cies, for several eigenstates of the three-fermion model with 
L = 15 and V = 10. Same plot in semilog scale is shown in 
the inset, indicating exponential localization. 



VI. DYNAMICS (TIME DEPENDENCE) 



It is instructive to study the localization phenomenon 
through explicit time evolution calculations. For the 
Bose-Hubbard model, some dynamical results have ap- 
peared in Ref . [l^ , In Figure [121 we present temporal 
dynamics results for the spinless-fcrmion model. 

For a three-fermion system in L = 10 sites, we present 
the evolution of the occupancies of individual sites, nj{t), 
after the system is started with the three fermions in 
consecutive sites. The left panels display the behavior 
when the fermions in the initial state are concentrated at 
one edge (the first three sites). There is no appreciable 
dynamics at the time scales shown, in marked contrast 
to the cases where the fermions start off at three other 
consecutive sites (middle panels and right panels). This 
difference is a dynamical demonstration of the edge lo- 
calization phenomenon. 
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FIG. 11: (Color online.) 1/L-dependence of max(nj) for (a) 
the eigenstates of the two-fermion chain, and (b) different 
eigenstates of the three-fermion chain. Here V = 10, and L 
goes from L = 9 to L = 21. Here S3 is the Hilbert space size. 



while the two exceptions are the edge states for which 
max{(nj)i,} has the constant value of 1/2. This demon- 
strates, once again, that all states in the two-fermion 
case are extended in space, while the three-fermion sys- 
tem possesses two localized states. 



Additional physics can be gleaned from the tempo- 
ral dynamics shown in Figure 1121 In the cases where 
we start from three non-edge consecutive sites (center 
and right panels), the first and last sites are never ex- 
cited {nj=i and Uj^iQ remain practically zero). This re- 
flects the spectral separation of the edge states, \e±) « 
(IlllOOOO...) ± |...0000111))/V2, from the rest of the 
breather band. Since the initial configuration is within 
the breather band, the dynamics is dominated by this 
band. Since the edge-localized subspace spanned by the 
basis {IlllOOOO...}, |... 0000111)} is separated from the 
rest of the breather band, the last/first sites are not ex- 
cited. 

Time scales. — Figure [T^ highlights dynamics at the 
scale of hundreds to thousands of t^^ units. The rea- 
son is that the dominant dynamics for our chosen initial 
conditions is that within the breather band; hence we 
are interested in coherent hopping of the three fermions. 
Such a cooperative hop event of three fermions, from one 
soliton (neighboring) configuration to the next, involves 
two intermediate states that are energetically ~ V away, 
in the (2-1-1) band. Thus the energy involved is of order 
t^/V^, so that the time scales are ^V^/t^. For V = 20, 
this leads to the above-mentioned time-scales. Perform- 
ing simulations at other values of V, we have seen that 
the relevant time scale indeed varies as ^V^. 

There is, of course, additional dynamics at other time 
scales. Tiny high-frequency wiggles can be seen in our 
data at scales of ^t~^, representing high-energy inter- 
band processes. Also, considering the left panel of Figure 
[T^ we note that the state 1 1110000...) is not itself an 
eigenstate; the true eigenstates are |e± ). Thus, there will 
be oscillations involving these two edge states, which is 
not visible here because the relevant time scales are much 
higher, and grow exponentially with the chain length. 
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FIG. 12: (Color online.) Dynamics of three fermions in ten sites, shown through the evolution of site occupancies rij. The 
three different initial conditions are shown on top. In each case, to avoid clutter we split the ten Uj into two groups of five. 
Thus the left top and left bottom describe two parts of the same evolving chain. Here V — 20 and t — 1. 



VII. CONCLUSIONS 

In this article, we have presented a straightforward and 
natural mechanism for localization in quantum lattice 
systems, namely the existence of open boundaries. Our 
edge-localization is a cooperative, as opposed to single- 
particle, phenomenon. We have provided perturbative 
arguments to explain the energy spectrum structure that 
lies at the heart of this localization phenomenon, thereby 
also explaining why at least three particles are required 
for the localization. In addition, we have showed the ap- 
pearance of richer edge configurations that appear for 
larger numbers of particles, and characterized the en- 
ergy spectrum and localization through a study of ex- 
plicit temporal dynamics. 

Energy scales. — At strong interactions, the bands 
are the most pronounced feature of the energy spectrum. 
The band energies are set by the interaction strength 
{U or V). One can then think of finer features of these 
bands, in various orders of the hopping t. There is always 
an energy shift at second order in i, since a basis state can 
be connected to itself through two hops. Whenever the 
degeneracy lifting requires higher than second order in t, 
we can have sub-band structures due to differing energy 
shifts at second order. In the two-boson or two-fermion 
case, this possibility does not exist. In the three-boson 
or three-fermion case, only the breather band splits at 
higher than second order, while the other two bands have 
degeneracy lifting at linear order. (The linear splitting 
is manifested in the shapes of the non-breather bands 
in Figures [T]-d and[5]-d.) Therefore in the three-particle 
systems only the breather band can have a subband; this 
subband turns out to have two members which are edge 



states. 

To rephrase, edge localization from the breather band 
can be seen as a result of the competition between (a) the 
distinct energy shifts ~ t'^/U recieved by the edge states, 
and (b) the degeneracy splitting for which the energy 
scale is ~ f^/W^^-^, for an n-boson system. For n > 3, 
analogous energy scale competitions can also play a role 
in other (multi-breather) bands where splitting occurs at 
larger order, leading to subbands and more complicated 
edge states. We have illustrated this for n = 4, in Section 

113 

Finally, we note that the two edge states separating 
out from the breather band are themselves very nearly 
degenerate at strong interactions. This degeneracy gets 
broken only at much higher order; the relevant energy 
scale is U{t/U)^ , with P ^ L for a chain with L sites. 

Interestingly, some of these energy scale issues have 
been probed dynamically in Section IVIl through the 
study of temporal evolution. 

Possible applications. — The edge localization has 
real-space effects on the dynamics, as revealed by our 
time evolution calculations (Section IVIl and Figure fT2)) . 
This raises the possibility of exploiting these effects for 
experimental quantum control of bosons or fermions in 
one-dimensional lattices. If one starts at a breather-band 
configuration that is not an edge state, the edge modes 
will not be excited, and conversely, fermions or bosons 
populated in an edge configuration remain stable in that 
configuration for a long time (Figure [T^ . One can con- 
ceivably use this effect, arising from fine structures in 
the energy spectrum, to manipulate and select sites for a 
few-particle lattice system in cold-atom or quantum wire 
experiments. 
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Open issues. — Our work opens up a host of open 
questions and issues, of which we mention a few. First, 
our calculations shown in Figure [4] indicate that for 
moderate-sized lattices, edge localization persists in some 
form at weak interactions. The fate of this weak- 
interaction localization for large lattice sizes remains an 
un-resolved question. 

Our work with two separate Hamiltonians results sug- 
gest that edge localization is a generic phenomena for 
quantum lattice models. Analogous phenomena can 
possibly be found in other one-dimensional itinerant 
fermionic and bosonic models. It is obviously of interest 
to investigate this phenomenon in various other quan- 
tum chain models, such as the one-dimensional fermionic 
Hubbard and extended Hubbard models. 

In the discrete nonlinear Schrodinger equation, one can 
find localized states at all lattice sites, both bulk and edge 
In contrast, for three or four bosons we have only 



found localization of all particles at the edge, even though 
some of our localized states have a finite width of more 
than one site. Can the Bose-Hubbard model support 
solitons some distance away from the edge, perhaps at 
larger boson numbers? 

Finally, since edge-localized classical breathers are 
known for two-dimensional lattices [2l|, there is the 
intriguing possibility of edge states in two-dimensional 
quantum lattices, which remains unexplored at present. 
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